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ABSTRACT 

We report here the results of the first Chandra X-Ray Observatory observations of the globular cluster M28 
(NGC 6626). We detect 46 X-ray sources of which 12 lie within one core radius of the center. We show that 
the apparently extended X-ray core emission seen with the ROSAT HRI is due to the superposition of multiple 
discrete sources for which we determine the X-ray luminosity function down to a limit of about 6 x 10-^° erg s '. 
We measure the radial distribution of the X-ray sources and fit it to a King profile finding a core radius of Cc.x ~ 11 
We measure for the first time the unconfused phase-averaged X-ray spectrum of the 3.05-ms pulsar B1821-24 and 
find it is best described by a power law with photon index E ~ 1.2. We find marginal evidence of an emission 
line centered at 3.3 keV in the pulsar spectrum, which could be interpreted as cyclotron emission from a corona 
above the pulsar’s polar cap if the the magnetic field is strongly different from a centered dipole. The unabsorbed 
pulsar flux in the 0.5-8.0 keV band is Ri 3.5 x 10“'^ergs s“* cm“^. We present spectral analyses of the 5 brightest 
unidentified sources. Based on the spectral parameters of the brightest of these sources, we suggest that it is a 
transiently accreting neutron star in a low-mass X-ray binary, in quiescence. Fitting its spectrum with a hydrogen 
neutron star atmosphere model yields the effective temperature = 90;!;[jj eV and the radius = 14.5;!;f g km. 
In addition to the resolved sources, we detect fainter, unresolved X-ray emission from the central core. Using 
the Chandra-derived positions, we also report on the result of searching archival Hubble Space Telescope data for 
possible optical counterparts. 

Subject headings: globular clustersigeneral — globular clustersdndividual (M28) — starsmeutron — x-ray:stars 
— binaries:general — pulsars:general — pulsarsundividual (PSR B1821-24) 


1. INTRODUCTION 

Since the Einstein era it has been clear that globular clusters 
contain various populations of X-ray sources of very different 
luminosities (Hertz & Grindlay 1983). The stronger sources 
(L^ w 10^®- 10^* ergs-*) were seen to exhibit X-ray bursts 
which led to their identification as low-mass X-ray binaries 
(LMXBs). The nature of the weaker sources, with Lx<3x 10^"* 
ergs"*, however was more open to discussion (e.g., Cool et al. 
1993; Johnston & Verbunt 1996). Although many weak X-ray 
sources were detected in globulars by ROSAT (Johnston & Ver¬ 
bunt 1996; Verbunt 2001), their identification has been difficult 
due to low photon statistics and strong source confusion in the 
crowded globular cluster fields, except for a few cases (Cool 
et al. 1995 and Grindlay et al. 1995). The application of the 
Chandra X-Ray Observatory (CXO) sub-arcsecond angular res¬ 
olution to the study of globular cluster X-ray sources leads one 
to anticipate progress in our understanding. 

Advances with the CXO have already included observations 
of NGC 6752 (Pooley et al. 2002a), 6440 (Pooley et al. 2002b), 
6397 (Grindlay et al. 2001) and 5139 (Rutledge et al. 2002a) 


which detected more low-luminosity X-ray sources than in all 
ROSAT observations of 55 globular clusters combined. Of par¬ 
ticular interest are the results obtained from CXO observations 
of 47 Tuc = NGC 104. Grindlay et al. (2001) reported the de¬ 
tection of 108 sources within a region corresponding to about 5 
times the 47 Tuc core radius. Sixteen of the soft/faint sources 
were found to be coincident with radio-detected millisecond 
pulsars (MSPs), and Grindlay et al. (2001,2002) concluded that 
more than 50 percent of all the unidentified sources in 47 Tuc 
are MSPs. This conclusion is in concert with theoretical esti¬ 
mates on the formation scenarios of short-period (binary) pul¬ 
sars in globular clusters (Rasio, Pfahl & Rappaport 2000). 

The globular cluster M28 = NGC 6626 lies close to the galac¬ 
tic plane (h = -5.58°) and close to the galactic center (/ = 7.8°) 
(Harris 1996). Distance estimates for M28 range from 5.1 kpc 
(Rees & Cudworth 1991) to 5.7 kpc (Harris 1996). In this pa¬ 
per, we use 5.5 kpc as a reference distance. 

M28 is a relatively compact cluster with a core radius of 
0(24, corresponding to ~ 0.4 pc, and a half-mass radius 
of 1(56, corresponding to ~2.6 pc (Harris 1996). The val¬ 
ues of these radii in parsec are smaller than those for the bet- 
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ter studied 47 Tuc. Thus, although M28’s central luminosity 
density, po = 10“*'^^ Lq pc“^, is comparable to that of 47 Tuc 
(10“* Lq pc“^ — Harris 1996), the rate of two-body encoun¬ 
ters in the core, A^enc c)C p^^rl, is a factor of 2 smaller, and thus 
fewer binaries are created and expected as dim X-ray sources 

Davidge, Cote & Harris (1996) noted that M28’s position in 
the plane near the galactic center and its relative compactness 
may indicate that it has been in the inner Galaxy for a long 
time. The authors suggest an age of ~16 Gyr. This age, how¬ 
ever, seems to be too high compared with the recent results by 
Salaris & Weiss (2002) and Testa et al. (2001) from which one 
estimates the age of M28 to be 11.4 -11.7 Gyr, consistent with 
M28’s moderately low metallicity [Fe/H]= -1.45. 

The absorbing column towards M28 is ~10 times larger than 
for 47 Tuc, with reddening E(B-V) = 0.43 (Harris 1996) corre¬ 
sponding to a hydrogen column density Nh - 2.4 X 10^' cm-2. 

The first millisecond pulsar discovered in a globular clus¬ 
ter was PSR B1821-24 in M28 (Lyne et al. 1987). This soli¬ 
tary MSP has a rotation period of P = 3.05 ms and period 
derivative of P = 1.61 x 10“'* s s“'. This value for the pe¬ 
riod derivative is sufficiently large that the expected correc¬ 
tion due to line-of-sight projection of acceleration in the clus¬ 
ter’s gravitational potential (Phinney 1993) is <10% (assum¬ 
ing a 10" projected offset from the cluster center and a clus¬ 
ter core mass ~ lO”* Mq; this assumed core mass is slightly 
higher than that derived from the values for central luminos¬ 
ity density and core radius given by Harris 1996). There¬ 
fore, the measured value reasonably accurately reflects the pul¬ 
sar’s intrinsic period derivative. The inferred pulsar param¬ 
eters make it the youngest (P/2P = 3.0 x 10^ yrs) and most 
powerful (E = 2.24 x 10*^/4g erg s“*) pulsar among all known 
MSPs. Here I 45 is the neutron star moment of inertia in units 
of 10"'^ g cm^. For the idealized magnetic dipole radiation 
model, the inferred perpendicular component of the magnetic 
dipole moment is /r sina = 3.2 x 10*^(/45PP)*/^ = 2.2 x lO^^/^g^ 
G cm*, where p is the magnetic dipole moment and a the an¬ 
gle between the rotation and magnetic dipole axes. Neglect¬ 
ing higher multipole contributions to the magnetic field and 
assuming the magnetic moment to be at the neutron star cen¬ 
ter, the inferred magnitude of the dipolar field at the magnetic 
pole is Bp = 4.5 x 10®/45^Pg*(sina)“* G, where Pg is the neu¬ 
tron star radius in units of 10® cm. Although these values for 
the dipole moment and polar dipole field are about 2 orders of 
magnitude smaller than what is inferred for ordinary field pul¬ 
sars in the galactic plane, they are the highest among all MSPs. 
The pulsar’s high rotational-energy loss made PSR B1821-24 
a prime candidate to be a rotationally powered, non-thermal 
X-ray source. This was confirmed by X-ray observations per¬ 
formed with different X-ray satellites, e.g., ROSAT (Danner et 
al. 1997 — cf. Verbunt 2001), ASCA (Saito et al. 1997) and 
RXTE (Rots et al. 1998), which detected Crab-like pulsations. 

While the detection of strongly pulsed emission would sug¬ 
gest a magnetospheric origin of the X-ray emission, a clear 
characterization of the PSR B1821-24 spectrum has been ham¬ 
pered so far by the crowding of sources in the region. Indeed, 
observation of M28 with the ROSAT HRI made it clear that 
all spectral data obtained from PSR B1821-24, especially in 
the “soft” bands <10 keV, suffer from spectral contamination 
from nearby sources. Stretching the angular resolution of the 
HRI to its limit, four X-ray sources within 2' of the cluster cen- 

^ Nenc OC (cf. Verbunt and Hut 1987) and Vdisp ~ (GM/rcf-^ cc p^^Tc, 


ter were discovered, including the two barely resolved sources 
RX J1824.5-2452E and RX J1824-2452P located within 15" 
of the center of M28, with the latter being identified from the 
timing as the counterpart of the millisecond pulsar B1821-42. 

Danner et al. (1997) put forth the suggestions that the ~ 
10" X 10" extended region of emission associated with RX 
J1824.5-2452E was due to either a synchrotron nebula powered 
by the pulsar or a number of faint LMXBs. The latter inter¬ 
pretation was favored as measurements with the ROSAT PSPC 
exhibited flux variations by as much as a factor of 3 between 
observations performed in 1991 and 1995 (see also Becker & 
Triimper 1999). Eurther evidence was provided by Verbunt 
(2001) who reanalyzed the ROSAT HRI data and resolved RX 
J1824.5-2452E into at least 2 sources. Time variations of the 
sources in M28 are also compatible with the result of Gotthelf 
& Kulkarni (1997) who discovered an unusually subluminous 
(Lx « 10“*LEdd) “Type I” X-ray burst from the direction of M28, 
thus leading one to expect the presence of one or more LMXBs. 

In this paper we report on the first deep X-ray observations 
of M28 using the ACIS-S detector aboard the CXO. We show 
that the X-ray core emission seen with the ROSAT HRI is dom¬ 
inated by a superposition of multiple discrete sources; we mea¬ 
sure the unconfused phase-averaged spectrum of the X-ray flux 
from the 3.05-ms PSR B1821-24; we establish the X-ray lumi¬ 
nosity function down to a limit of about 6 x 10*' erg s“' and 
measure, with high precision, the absolute positions of all the 
detected X-ray sources in M28 to facilitate identifications in 
other wavelength bands. Observations and data analysis are de¬ 
scribed in §2.1 through §2.3. In addition, we have used archival 
Hubble Space Telescope (HST) observations to search for po¬ 
tential optical counterparts of the sources detected by the CXO 
(§ 3 ). 

2. OBSERVATIONS AND DATA ANALYSIS 

M28 was observed three times for approximately equal ob¬ 
serving intervals of about 13 ksec between July and September 
2002 (Table 1). These observations were scheduled so as to be 
sensitive to time variability on time scales up to weeks. The 
observations were made using 3 of the CXO Advanced CCD 
Imaging Spectrometer (ACIS) CCDs (S2,3,4) in the faint timed 
exposure mode with a frame time of 3.241 s. Standard Chan¬ 
dra X-Ray Center (CXC) processing (v.6.8.0) has applied as¬ 
pect corrections and compensated for spacecraft dither. Level 2 
event lists were used in our analyses. Events in pulse invariant 
channels corresponding to « 0.2 to 8.0 keV were selected for 
the purpose of finding sources. Due to uncertainties in the low 
energy response, data in the range 0.5 to 8.0 keV were used for 
spectral analyses. Increased background corrupted a small por¬ 
tion of the third data set reducing the effective exposure time 
from 14.1 ksec to 11.4 ksec (Table 1) although no results were 
impacted by the increased background. 

The optical center of the cluster at 02000 = 18''24™32.®89 and 
^2000 = -24° 52' 11 "4 (Shawl and White 1986) was positioned 
1' off-axis to the nominal aim point on the back-illuminated 
CCD, ACIS-S3, in all 3 observations. A circular region with 
31 -radius, corresponding to twice the half-mass radius of M28, 
centered at the optical center was extracted from each data set 
for analysis. No correction for exposure was deemed necessary 
because the small region of interest lies far from the edges of 
the S3 chip. 

The X-ray position of PSR B1821-24 was measured sepa- 

where M oc por* is the mass of the core of the cluster (see also Verbunt 2002). 
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rately using the three data sets and the merged data. The results 
of these measurements are listed in Table 2. The set-averaged 
position is the same as that derived using the merged data set. 
The root-mean-square (rms) uncertainty in the pulsar position, 
based on the 3 pointings, is 0."042 in right ascension and 0"029 
in declination. The radio position and proper motion of the pul¬ 
sar, as measured by Rutledge et al. (2003), places the pulsar 
at the time of the observation only Aq = 0"083, = -0"042 

away from the best-estimated X-ray position. In what follows 
the observed X-ray positions of all sources have been adjusted 
to remove this offset. 

2.1. Image Analysis 

The central portion of the combined CXO image is shown 
in Figure 1. We used the same source finding techniques as 
described in Swartz et al. (2002) with the circular-gaussian 
approximation to the point spread function, and a mini¬ 
mum signal-to-noise (S/N) ratio of 2.6 resulting in much less 
than 1 accidental detection in the field. The correspond¬ 
ing background-subtracted point source detection limit is ~10 
counts. The source detection process was repeated using the 
CXC source detection tool wavdetect (Freeman et al. 2002) 
and yielded consistent results at the equivalent significance 
level. Forty-four sources were found using these detection al¬ 
gorithms. Close inspection of the image with the source posi¬ 
tions overlaid and of the source time variations showed that one 
source detected by the software was really two sources (num¬ 
bers 21 and 22 of Fig. 1), and that an additional source, # 24, is 
also present. 

Table 3 lists the 46 X-ray sources. The table gives the source 
positions, the associated uncertainty in these positions, the ra¬ 
dial distance of the source from the optical center, the signal- 
to-noise ratio, the aperture-corrected counting rates in various 
energy bands, the unabsorbed luminosity in the 0.5-8.0 keV 
energy range based on the spectroscopy discussed in § 2.2, and 
a variability designation according to the discussion in § 2.3. 

The positional uncertainty listed in column 4 of Table 3 is 
given by r = /N+Uo)^^'^ where a is the size of the circu¬ 

lar gaussian that approximately matches the PSF at the source 
location, N is the aperture-corrected number of source counts, 
and (To represents the systematic error. Uncertainties in the plate 
scale^® imply a systematic uncertainty of 0"13, and, given that 
the radio and CXO positions agree to 0."2, we feel that using 
0"2 is a reasonable and conservative estimate for (Jq. The fac¬ 
tor 1.51 is the radius that encloses 68% of the circular gaussian. 
The parameter a, varies from ~ 0."9 near the aimpoint to < 2f' 1 
near the edge of the 3(l-radius extraction region. 

The rates in the soft, medium, and hard bands listed in Ta¬ 
ble 3 are background subtracted and have (asymmetrical) 67% 
confidence uncertainties based on Poisson statistics (rather than 
the usual symmetrical Gaussian approximation). When the 
band-limited rates are positive, the uncertainties are symmet¬ 
rical in probability space (that is we set the lower and upper 
limits to the 67% confidence interval so that the true source rate 
is equally likely to fall on either side of the derived rate), but not 
in rate space. Otherwise, we set the rate to zero and calculate a 
67% upper limit. 

From Figure 1 we see that there are 12 point sources in the 
central region in the summed image. In addition, there re¬ 
mains some unresolved emission from this region of the cluster 


(§ 2.2.5). Pooley et al. (2002b) found similar diffuse emission 
in the central regions of the CXO image of the globular cluster 
NGC 6440. 

2.1.1. Radial Distribution 

The projected surface density, S{r), of detected X-ray sources 
was compared to a King profile. Sir) = 5o[l-l-(r/ro)^]“^-l-Co. 
The constant term Co was added to account for background 
sources. We estimated the number of background sources us¬ 
ing our observed flux limits of ~ 5.5 x 10“*® ergs cm“^ s“* in 
the 0.5-2.0 keV band and ~ 1.5 x 10“*® ergs cm“^ s“* in the 
2.0-10 keV band. These limits, together with the logA(> 5)- 
log5 distribution from the CXO Deep Field South (Rosati et 
al. 2002) gives an estimate of at least ~0.3 to 0.4 background 
sources per square arcmin in the field. (The true value may be 
higher because of the low galactic latitude of M28 relative to 
the deep field.) The best-fit value for Co is 0.36 ± 0.22 sources 
per arcmin^ indicating that there are 10 sources in the 3'- 
radius field not associated with the cluster. The other fit pa¬ 
rameters {So = 123 ±11 sources per arcmin^, ro = 23."8;!;g**9*, 
(3 = 3.5U[ ]) can be used to estimate the core radius, rc^x, and 
the typical mass, Mt, of the X-ray source population. We find 
a best fit core radius rc^x = lOf'O;!;* *. This is comparable to the 
distribution of optical light of the cluster: r„ = 14"4 (Harris 
1996). Following the derivation of Grindlay et al. (2002), the 
best-fit mass of the X-ray sources is Mt = l-87j;Q4g M©, as¬ 
suming the dominant visible stellar population has a mass of 
M* ~ 0.7 Mq. Although our range for Mt, estimated in this 
way, barely overlaps the range 1.1 -1.4 Mq deduced by Grind- 
lay et al. (2002) for 47 Tuc, additional effects due to uncertain¬ 
ties in cluster properties (position of cluster center, core radius, 
mass segregation, etc.) mean that the two results are in fact 
indistinguishable. 

2.2. Spectral Analysis 

Point-source counts and spectra were extracted from within 
radii listed in Table 4. Because the field is so crowded, back¬ 
ground was estimated using a region of 50" radius located 
in the south-western portion of the field. The background rate, 
used for all calculations, was 3.7 x 10“® counts s“* arcsec“^. 

Only 6 of the 46 detected sources have sufficient counts to 
warrant an individual spectral analysis. In descending order of 
the number of detected counts these are sources #26, #19, #4, 
#17, #25, and #28. Source #19 is the X-ray counterpart of the 
PSR B1821-24. The results of fitting various spectral models 
to the energy spectra of the brightest sources are presented in 
§ 2.2.1, § 2.2.2, and § 2.2.3. All spectral analyses used the 
CXC CALDB 2.8 calibration files (gain maps, quantum effi¬ 
ciency uniformity and effective area). The abundances of, and 
the cross sections in, TBABS (available in XSPEC v.11.2) by 
Wilms, Allen & McCray (2000) were used in calculating the 
impact of the interstellar absorption. All errors are extremes on 
the single interesting parameter 90% confidence contours. 

A correction** was applied when fitting models to the spec¬ 
tral data to account for the temporal decrease in low-energy 
sensitivity of the ACIS detectors due, presumably, to contami¬ 
nation buildup on the ACIS filters. The correction was based on 
the average time of the observations after launch of 1105 days. 

For the remaining 40 sources, all with fewer than 100 de¬ 
tected source counts, the spectra were combined to determine 


see http://asc.harvard.edu/cal/Hrma/optaxis/platescale/ 

* * available from http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/xanadu/xspec/models/acisabs.html 



4 


the mean spectral shape and total luminosity. Fits were made 
using both an absorbed power law d^V /dE oc E~^ (best-fit F = 
1.73 ±0.18, A^22 = A^/r/(1022cin-2) = 0.18 ±0.06, = 89.8 

for 86 degrees of freedom (dof)) and an absorbed thermal 
bremsstrahlung model {kT = 7.40±2.10keV, N 22 = 0.12 ± 0.04, 
X^ = 92.8). Using the power-law model parameters, the total 
(unabsorbed) flux from the 40 weak sources is 2.34 ±0.11 x 
10“^2 ergs cm“2 s“' in the 0.5-8.0 keV band, and the corre¬ 
sponding X-ray luminosity is 8.47 ±0.40 x 10^2 ergs s“'. This 
spectrum was then used to estimate the individual source lumi¬ 
nosities of the 40 faint sources listed in Table 3. See also the 
discussion in § 2.2.4. 

2.2.1. The Phase-averaged Spectrum of PSR 1821—24 

The spectrum of PSR B1821-24 was measured by extracting 
^1100 counts within a radius of l."72 centered on the pul¬ 
sar position. A background subtraction was performed, but its 
contribution (~ 2 counts) is negligible. Neither our source- 
detection algorithm nor wavdetect found any evidence for 
a spatial extent to the pulsar’s X-ray counterpart; thus, using 
the CXC model point spread function, 97% of all the events 
from PSR B1821-24 are within the selected region. The data 
were binned into 34 bins guaranteeing at least 30 counts per 
bin. Model spectra were then compared with the observed spec¬ 
trum. A power-law model was found to give a statistically 
adequate representation of the observed energy spectrum; the 
best-fit spectrum and residuals are shown in Figure 2. A ther¬ 
mal bremsstrahlung model resulted in a slightly better fit, but is 
not considered to be physically applicable. A blackbody model 
does not fit the data (x^ = 110 for 31 dof), as one could ex¬ 
pect based upon the hardness of the observed spectrum and the 
similarity of the sharp X-ray and radio pulse profiles (see, for 
example, Becker & Pavlov 2001, Becker & Aschenbach 2002). 

The best-fit power law yields N 22 = 0.16^QQg, F = 1.20^g [ 3 , 
and a normalization of 3.74^g 4 g x 10“^ photons cm“2 keV“^ 
at £ = 1 keV (xj = 0.89 for 31 dof). The column density is 
in fair agreement with what is deduced from the reddening to¬ 
wards M28. The unabsorbed energy flux in the 0.5-8.0 keV 
band is = 3.54);Qg5 x 10 “^^ ergs s“^cm“2, yielding an X-ray 
luminosity of = 1.28 ±0.02 x 10^2ergs s"^. This luminosity 
implies a rotational energy to X-ray energy conversion factor 
LxjE = 5.8 X 10““*. If transformed to the ROSAT band, this cor¬ 
responds to Lx = (3.4^.0) X 10^2ergs s“', and is similar to the 
luminosity inferred from the ROSAT data (Verbunt 2001). The 
photon index we found is compatible with F ~ 1.1 deduced for 
PSR B1821-24 from the observations at pulse maximum us¬ 
ing RXTE data (Kawai & Saito 1999). We note that we have 
ignored the possible effects of photon pileup, and this could 
artificially harden the spectral index. However, the degree of 
pileup here is sufficiently small (< 0.12 counts per frame), so 
that its effect on the spectrum is not significant. In fact, appli¬ 
cation of the Davis (2001) pileup model suggests the spectral 
index would be steeper by only 0.1 in the absence of pileup. 

The residuals in Figure 2 hint at a spectral feature or features 
at an energy slightly above 3 keV. Although the reduced num¬ 
ber of counts made it impossible to determine whether the fea¬ 
ture was present in the three separate observations, we note that 
some excess emission in the band between 3 and 4 keV is seen 
in all three data sets. By adding a gaussian “line” to the power 
law model, we found a line center at 3.3 keV with a gaussian 
width of 0.8 keV and a strength of « 6 x 10“® photons cm“2 
s“', which corresponds to a luminosity of « 1.1 x 10^' ergs“'. 


Adding the line changes xi to 0.63 (for 28 dof). The F-test 
indicates that the addition of this line component is statistically 
significant at 98% confidence, i.e., evidence for a broad spectral 
feature is marginal. 

If we assume that this feature is real, then it is interesting to 
speculate as to its origin. There are some atomic lines, from 
K and Ar, close to the 3.3 keV energy. However, K is not an 
abundant element, and it is hard to explain why the Ar lines are 
observed while no lines are seen from other, more abundant, 
elements. Therefore, we consider the more likely possibility 
that this is an electron cyclotron line, formed in a magnetic field 
B « 3 X 10 " G. Such a strong field can be explained by the pres¬ 
ence of either multipolar components or a strong off-centering 
of the magnetic dipole or both. 

For instance, if a dipole with the magnetic moment p = 
2 X 1022 G cm2 §1) is shifted along its axis such that it 
is 1.9 km beneath the surface, the magnetic field at the closest 
pole is 3 X 10 " G. The cyclotron line could be formed in an 
optically thin, hot corona above the pulsar’s polar cap, with a 
temperature kT ~ 10 keV such that the line’s Doppler width is 
smaller than observed while the excited Landau levels are popu¬ 
lated. The observed luminosity in the line can be provided by as 
few as Ag = 2.4 x 10^^ electrons, where T^ = T /(lO^K), 

Bii,5=B/(3 X 10"G). 

The corona is optically thin at the line center for an elec¬ 
tron column density tieH < 5 x 10'* s cm“2 (where H 

is the geometrical thickness of the corona) and a polar cap area 
Ape > 4.6 X 10® Tg 5 (polar cap radius Rpe > 12 m). 
The standard estimate of the polar cap radius, applied to the 
off-centered dipole, gives Rpe ~ r( 27 rr/cP )'/2 = 220 m, where 
r (= 1.9 km in our example) is the radial distance from the cen¬ 
ter of the magnetic dipole to the surface. At this value of Rpe, 

the optical thickness at the line center is 7" 0.003 7g^/^B72^. 

If the corona is comprised mainly of a proton-electron plasma, 
then its thickness can be estimated as H kT (mnsT^ ~ 40 Tg 
cm (at g ~ 2 X 10'^ cm s“2, typical for a neutron star), and a 
characteristic electron number density is ~ 5 x 10^* ^ 

cm“2. Thus, the cyclotron interpretation of the putative line 
looks quite plausible, and confirming the line with deeper ob¬ 
servations would provide strong evidence of local magnetic 
fields at the neutron star surface well above the “conventional” 
magnitudes inferred from the assumption of the centered dipole 
geometry. 

2.2.2. Source #26: An LMXB in quiescence? 

Of particular interest among the brightest X-ray sources de¬ 
tected in M28 is the luminous, soft source #26. We extracted 
data as per Table 4 and fit it with various spectral models. As 
expected from the softness of the count spectrum (Fig. 3), the 
power-law fit yields a photon index, F « 5.2, well in excess of 
those observed from known astrophysical sources with power- 
law spectra, and the hydrogen column density, N 22 ~ 0.68, that 
significantly exceeds the values expected from the M28’s red¬ 
dening (§1) and measured for the pulsar (§2.2.1). Therefore, 
we tried various models of thermal radiation. A blackbody 
fit gives the hydrogen column density N 22 ~ 0.13, tempera¬ 
ture A:7 bb ~ 0.26 keV and radius Rbb « 13 km, correspond¬ 
ing to the bolometric luminosity of Lbb ~ 1 x IO 22 erg s“'. 
Such values are typical for blackbody fits of LMXBs with tran¬ 
siently accreting neutron stars in quiescence (e.g., Rutledge et 
al. 2000). The relatively high temperatures of such old neutron 
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stars can be explained by heating of the neutron star crust dur¬ 
ing the repeated accretion outbursts (Brown, Bildsten, & Rut¬ 
ledge 1998). This heat provides an emergent thermal luminos¬ 
ity L = 8.7 X 10^^(M)_io erg s for the nuclear energy release 
of 1.45 MeV per accreted nucleon, where (M)_io is the time- 
averaged accretion rate in units of 1 O“^°M 0 yr“'. 

The radius Rbb obtained from the blackbody fit is much 
smaller than R ^ 10-15 km expected for a neutron star. A 
likely reason for this discrepancy is that the blackbody model 
does not provide an adequate description of thermal radiation 
from the neutron star surface. At the temperatures of interest 
the neutron star crust is covered by an atmosphere comprised 
of the accreted matter. Because of the gravitational sedimen¬ 
tation, the outermost layers of such an atmosphere, which de¬ 
termine the properties of the emergent radiation, are comprised 
of hydrogen, the lightest element present. Since the magnetic 
fields of neutron stars in LMXBs are expected to be relatively 
low, < 10® G, they should not affect the properties of X-ray 
emission, which allows one to use the nonmagnetic hydrogen 
atmosphere models (Rajagopal & Romani 1996; Zavlin, Pavlov 
& Shibanov 1996). We fit the observed spectrum with the nsa 
model'^ in XSPEC (v.l 1.2). In applying this model, we set the 
neutron star mass to I.TMq, leaving the radius of the emitting 
region and the surface temperature as free parameters. We ob¬ 
tained a statistically acceptable fit (xl, = 0.96, v = 44), with the 
hydrogen column density N 22 = 0.26 ± 0.04 (consistent with the 
expected value), the effective temperature kT^ = 0.09;!;Qg[ keV, 
a factor of 3 lower than 7 bb, and the radius = 14.5^3 g km, 
comparable with a typical neutron star radius. (The superscript 
°° means that the quantities are given as measured by a dis¬ 
tant observer; they are related to the quantities as measured at 
the neutron star surface as = grT, R°° = g~^R, 
where gr = {\-2GMiRc^yl^ is the gravitationalredshift factor. 
The values of and R°°, as inferred from the atmosphere 
model fits, are considerably less sensitive to the assumed value 
of neutron star mass.) 

The corresponding bolometric luminosity is = 1.9^g g x 
10^3 erg s '. Such a luminosity can be provided by a time- 
averaged accretion rate of (M) ~ 2 x 10 “^'Mq yr“'. Because 
all the fitting parameters we obtained for source #26 are typi¬ 
cal for quiescent radiation of other transiently accreting neutron 
stars in LMXBs (e.g., Rutledge et al. 2002b, and references 
therein), we conclude that such an interpretation is quite plau¬ 
sible. 

In addition to the thermal (photospheric) component, some 
X-ray transients in quiescence show a power-law high-energy 
tail (Rutledge et al. 2002b), apparently associated with a resid¬ 
ual low-rate accretion. We can see from Figure 3 that our ther¬ 
mal model somewhat underestimates the measured flux above 
2.5 keV. This is a consequence of a mild (0.16 counts per frame) 
amount of pileup. The pileup model suggested by Davis (2001), 
as implemented in XSPEC v. 11.2, does reproduce this tail, but 
only marginally changed the best-fit parameters. 

The quiescent emission of some of transient LMXBs shows 
appreciable variations of X-ray flux, with a time scale of a 
month (Rutledge et al. 2002b), and our source #26 also appears 
to be time variable (§2.3) in the three observations. We per¬ 


formed spectral fits to each of the data sets separately to check 
for spectral variation and found none, indicating that the time 
variability (§ 2.3) is not dominated by spectral variations. 

To test alternative interpretations of source #26, we also fit 
various model spectra of an optically thin thermal plasma in 
collisional equilibrium, applicable to stellar coronae and similar 
sources. Fitting the spectrum with the XSPEC model mekal 
for Z = 0.02Zq, we obtain a good fit (xj =0.88 for 44 dof) 
with N 22 ~ 0.33, kT w 0.6 keV, and emission measure EM 
w 2 X 10^^ cm“^, at d = 5.5 kpc. The corresponding luminos¬ 
ity, Lj:(0.5-8keV) w 1.2 x 10^^ erg s \ strongly exceeds the 
maximum luminosities of coronal emission observed from ei¬ 
ther single or multiple nondegenerate stars of any type avail¬ 
able in an old globular cluster. The spectrum is too hard, and 
the luminosity is too high, to interpret this emission as pro¬ 
duced by a non-magnetic CV (Warner 1995, and references 
therein). The inferred temperature is too low in comparison 
with the typical temperatures, ~ 30-40 keV, of the hard X- 
ray (bremsstrahlung) component observed in polars (magnetic 
CVs, in which rotation of the accreting magnetic white dwarf 
is synchronized with the orbital revolution), and, in addition, 
this component is usually much less luminous in polars. On the 
other hand, the spectrum of source #26 is too hard to be inter¬ 
preted as a soft X-ray component (20 - 40 eV blackbody plus 
cyclotron radiation) observed in many polars. Luminosities up 
to ~ 1 X 10^^ erg s“' (in the 2-10 keV band) have been ob¬ 
served in a number of intermediate polars (asynchronous mag¬ 
netic CVs with disk accretion). However, spectra of intermedi¬ 
ate polars are, as a rule, much harder (similar to those observed 
in polars) and strongly absorbed {N 22 ~ 10 ) by the accreting 
matter (Warner 1995). Therefore, we conclude that the inter¬ 
pretation of source #26 as a stellar corona or a CV looks hardly 
plausible, and most likely its X-ray emission emerges from the 
photosphere of a quiescent transient neutron star in an LMXB. 

2.2.3. Spectra of Four Other Bright Sources 

In addition to PSR B1821-24 and source #26, there are 4 
moderately bright sources with enough counts to attempt spec¬ 
tral fitting. Data from an extraction radius large enough to en¬ 
compass a significant fraction of the source counts were binned 
into spectral bins to maintain a minimum number of counts per 
bin, background subtracted, and fit to various spectral models. 
The extraction radii, total number of extracted counts, estimated 
number of background counts, minimum number of counts per 
spectral bin, and the number of spectral bins are listed in Ta¬ 
ble 4. To characterize the spectra of these sources, we fit each 
of them with three popular models of substantially different 
shapes: a black-body, a power-law, and an optically-thin ther¬ 
mal emission model (mekal) with Z = O.O2Z0. The best-fit 
parameters for these models, including uncertainties, are given 
in Table 5. 

The brightest of the four sources is source #4. Its spectrum 
(Fig. 4) is too hard to consider it as a qLMXB. The power-law 
fit of the spectrum, with E w 1 . 6 , might be interpreted as arising 
from magnetospheric emission from an MSP, with a luminos¬ 
ity Lx ~ 9 X 10^^ erg s \ in the 0.5-8 keV range. However, 
the hydrogen column density inferred from the power-law fit, 
N 22 ~ 0 . 86 , considerably exceeds those estimated from the in- 


These models are based on the work of Zavlin et al. (1996), with additional physics to account for comptonization effects (Pavlov, Shibanov & Zavlin 1991) 

This and similar models include thermal bremsstrahlung and line emission as components. Such models become equivalent to the optically thin thermal 
bremsstrahlung at high temperatures and/or low metalicities, while they strongly differ from the bremsstrahlung in the opposite case, when the main contribution 
to the soft X-ray range comes from the line emission. Therefore, there is no need to consider the thermal bremsstrahlung fits if more advanced models for optically 
thin plasma are available for fitting. 
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terstellar reddening and the power-law fit of the PSR B1821-24 
spectrum. The large Nh and a large distance, 2'.53 from the 
cluster’s center, together with the power-law slope typical for 
AGNs, hint that source #4 could be a background AGN (notice 
that the Galactic HI column density in this direction is about 
0.19 X 10^^ cm“^ — Dickey & Lockman 1990). The black- 
body fit gives N 22 « 0.16, consistent with that obtained for PSR 
B1821-24. This fit indicates that source #4 might be a ther¬ 
mally emitting MSP, although the blackbody temperature, « 1 
keV, is surprisingly high, and the blackbody radius, w 60 m, is 
much smaller than expected for a pulsar polar cap (a standard 
polar cap radius is Rpc = 1 AR^J^P_ 2 ^ km, assuming a centered 
dipole, where Rg is the neutron star radius in units of 10 km, and 
P-2 is the pulsar’s period in units of 10 ms). The inferred tem¬ 
perature would become lower by a factor of two, and the radius 
would increase by an order of magnitude if we assume that this 
emission emerges from a polar cap covered by a hydrogen or 
helium atmosphere (e.g., Zavlin & Pavlov 1998), but still a tem¬ 
perature of a few million kelvins looks too high, and a radius of 
a few hundred meters is somewhat too small, for a typical MSP. 
If the possible variability of this source (see Table 3) is con¬ 
firmed by future observations, the interpretation of this source 
as an MSP can be ruled out. The mekal fit gives a tempera¬ 
ture, kT ~ 30 keV, and a luminosity, L;t(0.5-8keV) ~ 8 x 10^^ 
erg s“\ too high to be interpreted as coronal emission of sin¬ 
gle or binary (BY Dra, RS CVn) nondegenerate stars. The high 
temperature'"* is typical of CVs, and the high absorption col¬ 
umn, N22 ~ 0.77, could be interpreted as additional absorption 
by the accreting matter, but the luminosity is somewhat higher 
than observed for most CVs (cf. Grindlay et al. 2001). Thus, 
the spectral fits suggest that source #4 is likely a background 
AGN, but they do not rule out the interpretation that it is a CV 
or an MSP. 

The spectrum of source #17 (Fig. 5) is even harder than that 
of source #4. It can be equally well fitted with a power-law 
model, with F w 1.3, and a mekal model, with kT ~ 40 keV 
and Lj:(0.5-8keV) ~ 5 x 10^^ erg s“'. The blackbody fit does 
not look acceptable, because of the high = 1.73 (the model 
underestimates the number of counts below 1 keV and above 
5 keV) and unrealistically low N22 < 0.03. The hardness of 
the spectrum is inconsistent with source #17 being a qLMXB, 
while its variability with a timescale of years (see §2.3) rules 
out a MSP interpretation. The luminosity of source #17 is lower 
than that of #4 so the argument against the CV interpretation is 
not so strong. On the other hand, source #17 is substantially 
closer to the cluster’s center, so the probability that it belongs 
to the cluster is higher. Therefore, a CV interpretation looks 
more plausible. 

Source #28 shows a softer spectrum (Fig. 6 ), in compari¬ 
son with sources #4 and #17, with possible absorption or an 
intrinsic turnover at softer energies. Because of the small num¬ 
ber of counts detected, we cannot distinguish between different 
fits statistically. Both the power-law fit (F « 3, N22 ~ 1-8) and 
mekal (kT k, 2 keV, N22 ~ 1-4) require an absorption column 
much higher than expected for a Galactic source in this direc¬ 
tion. On the other hand, the blackbody fit yields a lower (albeit 
rather uncertain) absorption, N22 ~ 0.3- 1.1. The blackbody 
temperature, kT^^ ^ 0.7 keV, and the radius, Rbb ~ 100 m, in¬ 
dicate that it might be thermal emission from an MSP polar 
cap. As we have discussed for source #4, a light-element atmo¬ 
sphere model would give a lower temperature, kT ~ 0.3 keV, 


and a larger radius, R ~ 1 km, which makes the MSP interpre¬ 
tation even more plausible. The observed spectrum of source 
#28 is harder than observed from qLMXBs, so we consider the 
qLMXB interpretation unlikely. 

Finally, the hard spectrum of source #25 (Fig. 7) strongly re¬ 
sembles that of source #17, although with much fewer counts. 
The blackbody model is unacceptable, while both the power- 
law and mekal yield reasonable fits. Similar to source #17, we 
consider source #25 as a plausible CV candidate. 

2.2.4. X-ray color-luminosity relation 

Figure 8 shows a plot of X-ray luminosity versus an X-ray 
“color”. Such X-ray “color-magnitude diagrams” (CMDs; see 
Grindlay et al. 2001) are particularly useful for studying source 
populations in clusters where a large dynamic range of source 
luminosities and types can be studied at a common distance. 
The source with the highest luminosity, and hardest spectrum, 
is the PSR B1821-24. The source with the second highest lu¬ 
minosity we consider to be a good candidate for a quiescent 
low-mass X-ray binary (qLMXB) based on its spectral prop¬ 
erties discussed in § 2.2.3. The other sources presumably are 
a mix of CVs (especially those with luminosities > 10^^ erg/s), 
RS CVns, main-sequence binaries, MSPs, and other (unknown) 
systems. 

Allowing for error bars on derived source colors and lumi¬ 
nosities, it is clear that the CMD is most useful for the classi¬ 
fication of brighter sources. However, even for faint sources, 
with poor statistics on each, a sufficiently large number of ob¬ 
jects in the CMD can define an approximate distribution of 
source types (e.g. the MSPs in 47 Tuc) and constrain the pos¬ 
sible source types of unidentified sources. In this observation, 
this has not been the case once the uncertainties are accounted 
for. Nevertheless, it is interesting to speculate that many of the 
soft, faint sources are lower luminosity MSPs as in 47 Tuc. 

2.2.5. Spectroscopy of the central unresolved emission 

As discussed in § 2.1, unresolved X-ray emission is present 
in the CXO data. This emission extends to roughly one core ra¬ 
dius from the center of M28. If we take the total counts detected 
within 15" and subtract off the known contribution from point 
sources (including the estimated counts from the full PSF) then 
we are left with an excess of about 541 counts. To extract this 
spectrum, counts near the point sources were removed. The re¬ 
sulting spectrum contained 568 counts in good agreement with 
our estimate for the excess. The background contribution to 
this total is ~100 counts. The spectrum was modeled using 
both a power-law model and a compound model consisting of a 
power-law with an additional optically-thin thermal emission¬ 
line (mekal) model. To account for the underabundance with 
respect to solar expected in the globular cluster, the abundance 
of metals in the mekal model were set to 2 % of their solar val¬ 
ues. The best-fitting model is the compound model (x^ = 30.4 
for 32 dof) though the mekal component is significant at only 
the ~ 2 a level (x^ = 37.3 for 34 dof for the power-law-only 
model). The best-fit photon index is F = 1 •79;!;o 32 , *he tem¬ 
perature of the emission-line component is kT = 0.18^q q* keV. 
The luminosity is 2.9 ;|;q 3 x 10^^ ergs s“' (6.0^['g x 10^^ ergs s“' 
after correction for absorption). 

Interestingly, the photon index is similar (but not identical) to 
that deduced from co-adding the 40 weakest resolved sources 


At such high temperatures the mekal model is essentially equivalent to the optically thin thermal bremsstrahlung. 
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(§2.2). This suggests a portion of the unresolved emission may 
be from point sources below the detection threshold but with 
similar spectral properties as those above the detection thresh¬ 
old. 

We find that the X-ray logA^(> 5)-log5 distribution of the 12 
sources within 15" of the cluster center is N(> S) ^ 525“® ^^. 
Assuming that this relationship extends to lower counting rates, 
we can estimate the number of photons that could come from 
sources below our threshold of 10 counts. This extrapolation 
predicts that ~ 200 counts or ~ 1.2 x 10 ^^ ergs s“' is con¬ 
tributed by sources below the detection threshold. Thus, un¬ 
resolved sources with spectral properties similar to the weaker 
resolved sources can account for roughly half the observed 
power-law component of the unresolved emission. 

We know that there are at least 4 distinct populations that 
could account for the unresolved emission: CVs, MSPs, 
BY Dra and RS CVn binaries, and isolated stellar coronae. 
If we assume that the unresolved emission, about 3 x 10^^ 
ergs s“^ (unabsorbed), is entirely due to stellar coronae and 
there are about 10^ stars in the volume that gives rise to the 
unresolved emission, then this implies that the average stellar 
corona radiates at ~ 3 x 10^^ ergs s“^. Since our Sun’s X-ray 
luminosity varies between 3 x 10^6 and 5 X 10^^ ergs s ' (Peres 
et al. 2000), we can say with confidence that the average star 
in M28 is less active than our Sun at its peak. This is not un¬ 
expected, since X-ray activity appears to correlate with rota¬ 
tion, and the old stars in M28 are likely slowly rotating. There¬ 
fore, we appeal to the usual suspects, faint CVs, MSPs, and RS 
CVn and BR Dra binaries, to account for the excess background 
emission. The X-ray luminosity functions for these classes are 
of course uncertain at the faint end. Assuming average lumi¬ 
nosities in the range 1 x 10 ^® to 2 X 10^° ergs s implies 50- 
1000 such sources. 

2.3. Time Variability 

We have used the three available ACIS observations to search 
for source variability on a time scale of weeks. The sepa¬ 
rate images are shown in Figure 9, while Figure 10 shows the 
variation in the counting rates. For all of the 46 sources, the 
number of counts detected in each of the 3 observations was 
computed by modeling the spatial distribution of events with 
a two-dimensional circular gaussian function with a width in¬ 
creasing according to the off-axis angle to approximately match 
the point spread function. Contributions from nearby sources 
were accounted for by using the same gaussian function out to 
a distance of 7 ct from the source considered. In each case, both 
the position and width of the gaussians were kept constant in 
the data fitting, with only the normalizations left as free param¬ 
eters. This approach gives a reasonable estimate for the number 
of counts from each source, even if the source was below the de¬ 
tection threshold and/or confused with another nearby source. 
As a check, the merged data was also examined and the results 
compared to the sums obtained from the 3 separate data sets. 
The observed counting rates agreed to better than 1.5%, indi¬ 
cating no obvious biases in our procedure. 

For each source, we then calculated the deviations of the 
number of counts measured in each observation with respect 
to a constant flux distribution. By adding the 138 measure¬ 
ments (46 X 3), this calculation yielded a of 379.7 for 92 dof, 
clearly showing that some of the sources actually varied. The 
second observation of source #10 contributed the most (44.7) to 
X^, and so we designated that source as “variable”. The three 


observations of this source were then removed from the sample. 
We then repeated this process. In this way we determined that 
12 of the 46 sources exhibit evidence for some form of time 
variability. Setting aside these 12 sources, was 128.4 for 
102 measurements and 68 degrees of freedom. The largest con¬ 
tribution now (#28, 2nd measurement) contributes only 5.95 to 
X^ (2.44a deviation). Although it is likely that this source also 
varies (as do a few others) the level is low enough to make char¬ 
acterizing the variability difficult, and thus we did not ascribe 
any designation of variability to these sources. 

To characterize the variability further, if a source in one ob¬ 
servation is more than a factor of 2 above the average for the 
other two measurements, we say that it brightened ("b" in Ta¬ 
ble 3). Likewise, if one observation is more than a factor of 2 
below the average for the other two measurements, we say that 
it dimmed ("d" in Table 3). If the source neither “brightened” 
nor “dimmed” according to these definitions, but was still iden¬ 
tified as varying we denote it with a "v". 

The source identified as a possible qLMXB (#26) was 
flagged as variable. Since the emission is dominated by a very 
slowly cooling black body (§2.2.2), variability is not expected. 
The variability designation is a consequence of the third obser¬ 
vation being 13 percent (4a) higher than the first two. There 
are a number of possible systematic effects that might produce 
such an effect. We verified the integration time for all three 
observations by counting the number of ACIS frames. We 
also searched for missed bad pixels and columns by checking 
which physical pixels were included in the extraction region, 
and found none. Furthermore, the small change in the space¬ 
craft roll (12 degrees) between observations combined with the 
spacecraft dither, resulted in a common set of pixels for all three 
measurements. It is worth noting that cosmic ray tracks, the ef¬ 
fects of which are discarded by the ACIS flight software, can 
produce up to 10 percent variations in sensitivity, but of the 
front-illuminated (FI) chips. Our observations, however, were 
performed with a back-illuminated CCD. These are much thin¬ 
ner than the FI chips and consequently these tracks are about 
10 times smaller, with a corresponding smaller impact on sen¬ 
sitivity. In summary, we were unable to find any instrumental 
effect that could account for the observed variability for source 
#26. On the other hand, the small amplitude of the variability 
suggests caution in applying the designation. 

In addition to the millisecond pulsar, we were also able to 
compare counting rates with one other ROSAT observation — 
that of the source with the fourth largest number of detected 
counts, #17. This source is located outside of the core radius at 
42" from the cluster’s optical center. Although this source ap¬ 
pears constant during our CXO observations, ROSAT observa¬ 
tions indicate that the source is variable on time scales of years. 
Specifically, this source was not detected in either the ROSAT 
PSPC observation of March 1991 nor the ROSAT HRI obser¬ 
vation of September 1995. The source was, however, detected 
with ROSAT HRI in September 1996. 

To summarize, in the CXO data we found 6 sources that 
brightened, 3 that were variable, and 2 that dimmed. In ad¬ 
dition, we found one variable source from the ROSAT data. 
Seven of these 13 sources are within one core radius of the cen¬ 
ter of the globular cluster. Of these seven, the three that bright¬ 
ened are most probably main sequence stars, whereas the two 
that were designated variable and the two that dimmed are most 
likely CVs (or perhaps qLMXBs). Finally, all the sources that 
appear to vary on shorter time scales (within one of the obser- 



vations) had already been identified as varying by the technique 
described above. 

3. OPTICAL OBSERVATIONS 

We have performed a search for potential optical counterparts 
of the X-ray sources listed in Table 3 using data obtained with 
the HST WFPC2 and available in the public HST data archive. 
The observations of the M28 field were taken with the “V-band” 
filter F555W (A = 5500 A; AA = 1200A) and the “I-band” fil¬ 
ter F814W (A = 7995 A; AA = 1292 A) on 1997 September 12 
(Testa et al. 2001). To allow for a better cosmic ray filtering, the 
observations were split into a sequence of eight 140 s exposures 
in the F555W filter, and three 180 s plus six 160 s exposures in 
the F814W filter. Three short exposures of 2.6 s each were 
acquired in both filters to obtain unsaturated images of bright 
cluster stars. The total integration time was 1130 s and 1510 s 
in the F555W and F814W filters, respectively. The same data 
have been used by Golden et al. (2001) to search for the optical 
counterpart of the MSP B1821-24. 

Data reduction and photometric calibration were performed 
through the HST WFPC2 pipeline. For each filter, single expo¬ 
sures were combined using a cosmic ray filter algorithm. The 
final images were then registered on each other. Automatic 
object extraction and photometry was run by using the RO- 
MAFOT package (Buonanno and lannicola 1989). The source 
lists derived for each passband were finally matched to produce 
the color catalog. Conversion from pixel to sky coordinates 
was computed using the task metric which also applies the cor¬ 
rection for the WFPC2 geometrical distortions (see Testa et al. 
2001 for further details on the data reduction and analysis). 

The final catalog, consisting of a total of 33972 entries, with 
coordinates and magnitudes in the F555W passband and the 
F555W-F814W color, has been used as a reference for the op¬ 
tical identification. As a first step, the optical catalog has been 
cross-correlated with our list of X-ray sources. Since the ACIS 
astrometry has been boresighted using the pulsar’s radio coor¬ 
dinates as a reference (§2.1), the cross-correlation radius ac¬ 
counts only for the statistical error of the X-ray position and for 
the uncertainty of the HST astrometry. The latter is ascribed 
to the intrinsic error on the absolute coordinates of the GSCl.l 
guide stars used to point the HST and compute the astrometric 
solution. Typical uncertainties are of the order of 1"0 (see, e.g. 
Biretta et al. 2002) 

Twenty-two of the X-ray sources were in the HST field- 
of-view, and each of them has potential WFPC2 counterparts. 
Figure 11 shows the location of all the 376 matched WFPC2 
sources on the color-magnitude diagram (CMD) derived from 
all the sources detected in the WFPC2 field of view. The mea¬ 
sured magnitudes have been corrected for the interstellar red¬ 
dening assuming the same color excess E(B-V) = 0.43 (Harris 
et al. 1996) for all the sources. Most of the candidate coun¬ 
terparts lie on the globular cluster main sequence, with only 
a few of them possibly associated with evolved stellar popu¬ 
lations. We note that using any other average value of the red¬ 
dening [e.g., E(y-/) = 0.52, derived directly from our datasets], 
merely shifts the whole CMD. However, a more accurate anal¬ 
ysis of the stellar populations may require accounting for pos¬ 
sible differential reddening along the line of sight (Testa et 
al. 2001). This could narrow the distribution. 

Additional F555W (340 s) and F814W (340 s) WFPC2 ob¬ 
servations of the M28 field, taken on 1997 August 8, have been 

*5 www.stecf.org 


used to search for variability among the potential WFPC2 coun¬ 
terparts. The data have been retrieved from the ST-ECF pub¬ 
lic archive'^ after on-the-fly recalibration with the best avail¬ 
able reference files. The two pointings are centered very close 
to each other but with a slight relative rotation angle (~ 5° ). 
Image co-addition, object detection and position measurements 
were performed consistent with the previous analyses. To avoid 
systematic effects due to the difference in the default astromet¬ 
ric solution between the two datasets, which has been measured 
to be of the order of l."0, the object catalogs derived from the 
August 8 and September 12 observations were matched in the 
pixel space after registering the two sets of coordinates through 
a linear transformation. The overall dispersion of the radial co¬ 
ordinate residuals after the transformation turned out to be 0.17 
WEC pixels (0."017). Eor this reason, only objects with radial 
coordinate residuals smaller than 0.5 WEC pixels, i.e. 3 cr of 
the dispersion of the residuals, were considered as matched in 
the two datasets and examined for variability. Possible spurious 
matches were checked manually and filtered out. 

Although a number of candidates show brightness variation 
larger than 0.5 magnitudes in at least one of the filters, their na¬ 
ture cannot be assessed with high confidence. Some of them are 
indeed close to the detection limit, hence with larger photom¬ 
etry errors, while others are detected too close to bright stars 
to obtain clean measurements even using PSP subtraction al- 
ghoritms, and a few are located in very crowded patches. 

We note that some objects detected in the first dataset are ab¬ 
sent in the second. All these cases have been checked carefully 
to find out whether the lack of matches was due to intrinsic ob¬ 
ject variability. However, we found that the missing matches 
can be explained either by the fact that in the second dataset 
some objects fall at the chip edge or in the overscan region, or 
they fall out of the field of view or, as the second dataset is 
shallower, they are simply below the detection limit. 

Although we have found a number of candidate counterparts 
to the X-ray sources, no definite conclusions can yet be drawn 
from our search. In the crowded globular cluster field, the dom¬ 
inant uncertainty in the HST astrometry represents the real bot¬ 
tleneck for obtaining optical identifications. In contrast, e.g., to 
the experience of Pooley et al. (2002a) with NGC 6752, we 
found very few blue candidates which can be considered as 
likely counterparts and used to Foresight the HST astrometry. 
At the same time, no GSC2.1 or USNO-A2.0 stars are present 
in the narrow WPPC2 field of view to use them as a reference to 
recompute the image astrometric calibration. The only way to 
improve the WPPC2 astrometry of our datasets is by upgrading 
the coordinates and the positional accuracy of the guide stars 
used for the telescope pointing and recomputing the astromet¬ 
ric solution in the HST focal plane. This, together with a larger 
spectral coverage, will definitely give us a better chance to ob¬ 
tain firm identifications. 

4. SUMMARY 

We have analyzed observations of the globular cluster M28 
taken with the ACIS-S3 instrument aboard the CXO. Porty-six 
X-ray sources were detected within 3.'1 of the optical center of 
the cluster down to a limiting (absorption-corrected) luminos¬ 
ity of ~ 6.8 X 10^*' ergs s“' in the 0.5 - 8.0 keV range. Many 
of these sources are concentrated near the center of the clus¬ 
ter. Their radial distribution can be described by a King profile 
with an X-ray core radius comparable to that of optical light 
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but with a steep power-law index suggesting a rather high X- 
ray source population mass, ~ 1.9 Mq. The X-ray source 
distribution flattens at larger radii consistent with a population 
of background sources. 

Among the brightest sources in M28 is the millisecond pulsar 
PSR B1821-24. We find the phase-averaged spectrum of PSR 
B1821-24 to be best represented by a F = 1.20 power law ra¬ 
diating at Lj: = 1.3 X 10^^ ergs s“'. The luminosity is consistent 
with a steady luminosity since the time of ROSAT observations. 

An intriguing spectral feature, albeit of marginal statistical 
significance, observed at ~3 keV in PSR B1821-24 might be 
an electron cyclotron line. If so, the line requires a magnetic 
field of order 100 times the strength inferred from P and P sug¬ 
gesting the local magnetic field at the surface of the neutron star 
is well above the conventional values obtained by assuming a 
centered dipole geometry. This may be due to multipolar com¬ 
ponents, to a strong off-centering of the magnetic dipole, or to 
a combination of the two. Recently, Gil & Melikidze (2002) 
argued on both observational and theoretical grounds that such 
strong deviations from the dipole magnetic field should exist at 
or near pulsar polar caps. Geppert & Rheinhardt (2002) showed 
that it is possible to create strong but small-scale poloidal field 
structures at the neutron star surface via a Hall-instability from 
subsurface toroidal field components. 

A second bright source, closer to the center of M28, was also 
studied in detail. The spectrum of this source (#26) is notably 
soft and thermal (Fig. 3) as is typical of transiently-accreting 
neutron stars in quiescence. Such objects have a hydrogen-rich 
atmosphere comprised of matter accumulated and heated dur¬ 
ing previous accretion episodes. Non-magnetic hydrogen at¬ 
mosphere models provide the best fit to the spectrum of this 
source and thereby support this interpretation. The bolomet- 
ric luminosity corresponding to the best fit model parameters, 
^ 1.9 X 10^^ ergs s“', can be maintained by a time-averaged 
accretion rate of ~ 2 x 10“" M© yr“'. While the flux from this 


source varies among the 3 CXO observations, there is no evi¬ 
dence for spectral variability between our observations taken at 
roughly 1 month intervals. 

While CXO resolves many of the X-ray sources in M28, 
there remains some diffuse emission distributed over ~1 core 
radius. This emission is only 22% of the total X-ray luminosity 
in the central region. Nevertheless, simple extrapolation of the 
observed logA - log5 relation to below our detection threshold 
cannot account for more than about one-half of the emission. 
Another population of weak sources must account for the re¬ 
mainder. Possibilities include BY Dra and (weak) RS CVn sys¬ 
tems, other millisecond pulsars, isolated low-mass stars, and 
CVs. 

By scheduling our 3 observations of M28 over an ~2 month 
period, we were able to observe variability in the resolved X-ray 
source population on a time scale of order weeks. Twelve of the 
sources, or 26%, are seen to vary over the course of the obser¬ 
vations including the quiescent low-mass X-ray binary. Some 
sources exhibited a much higher flux in 1 of the 3 observations, 
and may be associated with stellar coronae. Others display the 
opposite effect, being low in one observation or otherwise var¬ 
ied and may be associated with CVs and qLMXBs. 

Finally, the benefit of an accurate radio position for PSR 
B1821-24 has allowed us to constrain the X-ray positions of the 
M28 sources to sub-pixel accuracy. Comparison to HST images 
suffers, however, from the intrinsic error in the absolute coor¬ 
dinates of the GSCl.l guide stars used to compute the WFPC2 
astrometric solution thus preventing any conclusive identifica¬ 
tions. 

Those of us at the Marshall Space Flight Center acknowl¬ 
edge support from the Chandra Project. G.G.P. acknowledges 
support from NASA grant NAG5-10865. V.T. acknowledges 
support from MIUR under grant COFIN 2001: "Modeling old 
Millisecond Pulsars in the Galaxy and in Globular Clusters" 
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Fig. 1. — CXO ACIS-S3 image of the central region of M28. Twelve X-ray sources are detected within the 0f24 core-radius, indicated by a dashed circle. The 
optical center of the cluster is indicated by a cross. The X-ray counterpart of the millisecond pulsar B1821-24 is source #19. 




normalized co 




)ulsar PSR B1821—24 fit to an absorbed power-law model {upper panel) and contribution to the fit statistic 
significant feature near 3 keV. 



channel energy (keV) 


Fig. 3.— Energy spectrum of source #26. The spectrum is modeled with a nonmagnetic neutron star H atmosphere model. 
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channel energy (keV) 


ce #25, fit with an absorbed power-law model. 
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Fig. 8.— X-ray luminosity (0.5-8.0 keV), in units of 10^^ erg s“b vs. X-ray color for the sources listed in Table 3. The X-ray color (softness) is here defined as 
the rate in the 0.2-2.0 keV band divided by the rate in the 2.0-8.0 keV band. The open circle marks the softness for the average spectrum of the 40 faintest sources 
at their average luminosity of 2.25 x 10^^ erg s“^. Error bars are displayed for selected sources. For the dimmest sources. X-ray color is, practically speaking, 
unknown. The PSR B1821-24 is the most luminous source in our sample. The source at the top right of the diagram, #26, is possibly a qLMXB. 
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RIGHT ASCENSION (2000) 


Fig. 9.— ACIS-S3 image of the central region of M28 at three separate epochs (A, B, C) as per Table 1. 
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10 20 30 40 

Source number 


Fig. 10.— The count rate for all sources for each of the three observations. The x-axis is the row number listed in Table 3. For each source, three rates are plotted: 
a circle for the first observation, an “x” for the second, and a triangle for the third. Any missing points were below the minimum rate plotted. The statistical error 
bar is also plotted, but for the high count rate sources it is difficult to see. The variability indicators (b, v and d) are discussed in § 2.3. 
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Fig. 11.— Color-magnitude diagram for all the sources detected in the WFPC2 field of view (light grey points). The magnitudes are corrected for the interstellai* 
extinction using the same E(fi-F) = 0.43 for all the sources. The 376 HST candidate counterparts matched to the X-ray sources listed in Table 3 ai‘e indicated by 
grey filled circles. 
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Table 1 

CXO Observations 


Date 

ObsID 

Exposure 

2002 July 4 

2684 

12.7 ksec 

2002 Aug 8 

2685 

13.5 ksec 

2002 Sep 9 

2683 

11.4 ksec 


Table 2 

PSR B1821-24 Positions (J2000) 




Date 

Position 


2002 July 4 

18 

24 

32.015 

-24 

52 

10.81 

2002 Aug 8 

18 

24 

32.016 

-24 

52 

10.76 

2002 Sep 9 

18 

24 

32.009 

-24 

52 

10.83 

average 

18 

24 

32.013 

-24 

52 

10.80 

rms (arcsec) 


0 .( 

342 

0.029 

merged data 

18 

24 

32.013 

-24 

52 

10.80 

radio (8/02) 

18 

24 

32.008 

-24 

52 

10.76 
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Table 3 

M28 Discrete X-Ray Sources 



R.A. 

Dec. 


/ 

S/N 

Soft" 

Medium‘s 

Hard" 

Lx^ 

Variability^ 


(J2000) 

(J2000) 

(") 

(") 


(10“’ s“') 

(10“ s“') 

(10“ s-‘) 

(10”ergs“‘) 


1 

18 24 20.531 

-24 51 33.04 

0.433 

172 

4.09 

15.7!';" 


10 1+'^ 

17.0 


2 

18 24 20.619 

-24 51 27.17 

0.431 

172 

4.16 

15.8!’:J 

oc o+lO.O 

13.3+“ 

17.8 

b 

3 

18 24 22.575 

-24 52 05.69 

0.361 

140 

5.02 

35.5!|“ 

36.3!l“ 

663.4““ 

17.3!?| 

25.1 


4 

18 24 22.684 

-24 51 02.65 

0.303 

154 

21.12 

47.9!|‘;’ 

'jcfi c445.8 

756.5_46.2 

578.9 

V 

5 

18 24 22.831 

-24 52 46.06 

0.469 

141 

2.60 

13.2^^ 

Q 0+4.8 
^•^-10.6 
'3Z q+10.9 

36.8 116 

7 +8.^ 
^•-2.0 

6.5 


6 

18 24 24.237 

-24 51 04.14 

0.421 

135 

3.75 

14.9!“ 

20.l!|] 

1 ^0+8.8 

15.3 


7 

18 24 24.603 

-24 53 02.10 

0.392 

123 

3.14 

17 '^+7-6 

7 "7+11.1 

9.9 


8 

18 24 25.175 

-24 54 06.67 

0.327 

155 

5.75 

54 

55 5+*^-^ 

27 

33.4 


9 * 

18 24 25.189 

-24 52 13.60 

0.408 

104 

2.82 

044.2 

21 3‘'’^'^ 

044.2 

31.3^ 

36.5!|“ 

9.3 


10 

18 24 25.658 

-24 50 34.14 

0.407 

138 

4.58 

12.9“;’ 

20.4 

b 

11 

18 24 28.449 

-24 50 33.59 

0.423 

114 

3.72 

11 4-^^ 5 

OO Q+l5.i 

11.8“’ 

14.6 


12* 

18 24 28.580 

-24 53 09.77 

0.359 

82 

3.11 

0+3.2 

17.1+“ 

16.9+^j 

9.7 


13 * 

18 24 28.727 

-24 51 24.56 

0.405 

73 

3.23 

lO.t] 

13.4!“ 

10.2“;] 

10.2 


14* 

18 24 30.155 

-24 51 49.81 

0.323 

42 

5.76 

18.9!“ 

68.7!i“ 

48 7+^^-^ 

34.8 


15 * 

18 24 30.770 

-24 52 33.19 

0.367 

36 

2.91 

1 'X 7+6.2 
'-6.9 

14.1!“ 

0+3.2 

8.3 


16* 

18 24 30.946 

-24 52 13.84 

0.323 

26 

5.12 

^ 0+8.1 

8 

53.5!|“ 

07 0+9.5 

27.5 


17 * 

18 24 31.063 

-24 52 45.20 

0.299 

41 

18.77 

254.3!g;* 

537.6““ 

514.9!]^ 

439.1 

v" 

18 * 

18 24 31.591 

-24 52 17.49 

0.377 

18 

2.88 

5.3!^'’ 

14.5!“ 

11 3^^'^ 

7.8 


19* 

18 24 32.008 

-24 52 10.76 

0.298 

11 

30.09 

4+38.4 

1369+g 

1405+] 

96.0 


20 

18 24 32.213 

-24 53 51.58 

0.331 

100 

3.69 

21.5+^-^ 

24.5!“ 

13.4+^^ 

15.6 


21 * 

18 24 32.272 

-24 52 09.46 

0.323 

8 

4.85 

29.5+^-^ 

oo /c+14.7 

77.6_,6, 

47 «:+ll-5 
4/.D_,7 , 

27.6 


22 * 

18 24 32.345 

-24 52 08.02 

0.316 

8 

5.30 

35.4!|“ 

1 1 Z 4+18.0 
il6.4_i8 5 

£Q 4+i3.8 
68.4 _i4 4 

37.2 

d 

23 * 

18 24 32.514 

-24 52 07.66 

0.312 

6 

6.73 

44.3!|“ 

VD.O_i69 

137.3““ 

64.5 


24* 

18 24 32.631 

-24 52 21.70 

0.354 

10 

3.35 

1 ^ 7 + 6.2 
'-6.9 

22.7*11 

13.8“;] 

12.0 

b 

25 * 

18 24 32.689 

-24 52 23.54 

0.304 

12 

9.65 

54 +‘2-1 
•^^■- 12.6 

129.l!l“ 

172 

^ '"^-- 22.0 

129.0 


26* 

18 24 32.821 

-24 52 08.26 

0.298 

3 

36.93 

2213.8+™, 

2477. l!g’ 

211 . 7 ““ 

534.0 

V 

27 * 

18 24 33.026 

-24 52 01.72 

0.353 

9 

3.46 

8-7!5i‘3 

l^--6.7 

24.!“ 

12.7 

b 

28 

18 24 33.026 

-24 50 52.86 

0.306 

78 

13.72 

42.9!1?;5 

£1-1 ^+13.7 

O-TQ +27.6 

298.1+]]j 

207.9 


29* 

18 24 33.070 

-24 52 10.45 

0.308 

2 

7.90 

124.4+!|j 

q+13.7 

67.8_,4-, 

65.8 

V 

30* 

18 24 33.429 

-24 52 15.49 

0.313 

8 

6.50 

34.6!?„*3 

76.!|« 

C7 0+12.0 

42.8 

b 

31 

18 24 33.539 

-24 52 33.16 

0.316 

23 

5.62 

13.67J 

61.9!|“ 

0^ 0+9.8 

3°-3-10.3 

/r+13.7 

31.5 


32* 

18 24 33.634 

-24 52 12.12 

0.312 

10 

6.74 

17.2+“ 

103.4!l“ 

47.3 

d 

33 » 

18 24 33.759 

-24 52 11.08 

0.321 

11 

4.87 

OQ 1 +10.3 

71.1!|“ 

20 . 4 +“ 

29.4 


34* 

18 24 33.861 

-24 51 11.99 

0.391 

60 

3.59 

7.4!{o''2 

24.5!“ 

7 5‘^-^ 

11.8 


35 » 

18 24 34.469 

-24 53 12.99 

0.335 

65 

3.67 

26.9!“ 

2i.7!J:J 

0+3. 

^-0.0 

13.4 


36 

18 24 34.974 

-24 54 57.88 

0.368 

168 

3.39 

18.5!“ 

28.7“ 

12 . 2 “] 

11.3 


37 

18 24 35.524 

-24 52 18.34 

0.321 

36 

5.45 

14.3+t'J 

«--2 0 + 12.1 

35 2^^'^ 

29.5 


38 

18 24 36.578 

-24 50 16.04 

0.340 

125 

7.58 

52.9!|“ 

74 -1+14.2 

64.8+“] 

56.5 


39 

18 24 37.326 

-24 51 57.50 

0.363 

61 

3.46 

25.9!“ 

16 4+'^-^ 

7 7+i0.8 
Z./_2 7 

13.2 


40 

18 24 37.831 

-24 51 44.58 

0.401 

72 

2.93 

13.1+1, 

17.5!i^ 

0 £+1 1.1 
2-0_2,6 

9.1 

b 

41 

18 24 37.897 

-24 49 28.45 

0.632 

176 

2.75 

4, c+10.7 

Z.D_2_5 

19.4^ 

^•^-2.5 

9 . 1 “;] 

8.6 


42 

18 24 39.061 

-24 51 08.26 

0.434 

105 

3.15 

4 7+18- 
'^.1 

29 . 1 !“ 

1 2^''^ 
^■•^-1.3 

10.3 


43 

18 24 39.567 

-24 50 30.18 

0.561 

136 

2.63 

12 5“^'^ 

10.9!|S 

6.8 


44 

18 24 40.336 

-24 50 35.31 

0.563 

139 

2.64 

12.5+^’ 

4 7+16.3 

6.8 

b 

45 

18 24 41.303 

-24 54 16.97 

0.419 

170 

2.63 

044.1 

25 4+^-'* 

7 7+11.4 

8.1 


46 

18 24 42.694 

-24 52 47.48 

0.433 

138 

2.81 

10 1+^"^ 

1 1 /i+6-0 

13.6_,, 

4 7+ik2 

4./_4 7 

8.1 



Sources indicated with a * are in the field of view of the HST observations discussed in § 3 
Positional uncertainty radius in arcsec (see text). 

* Distance from the nominal optical center of the cluster. 

^ Detected counting rate, corrected for the psf, in the 0.2-1.0 keV band. 

^ Detected counting rate, corrected for the psf, in the 1.0-2.0 keV band. 

^ Detected counting rate, corrected for the psf, in the 2.0-8.0 keV band. 

^ X-ray luminosity in the 0.5-8.0 keV band assuming a distance of 5.5 kpc and Nh of 0.18 X 10^^ cm~^. The luminosities of the six brightest sources 
are based on a canonical power-law. More accurate luminosities for the six brightest sources are presented in the text. 

^ V indicates variability, f a flare and d a dip as discussed in § 2.3. 

^ based on comparison with ROSAT observations. 
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Table 4 

Extraction and spectral fitting parameters 


source # 

’ ext 

fi 


Ni 

N‘^. 

mm 

nL 

19 

l."72 

97.5 

1119 

2 

30 

34 

26 

l."72 

97.5 

1669 

2 

30 

47 

4 

8 ."9 

100 

540 

49 

15 

33 

17 

7."5 

100 

527 

28 

20 

24 

28 

14."8 

95 

300 

115 

15 

18 

25 

l."35 

95 

127 

1 

10 

11 



Note. — 


“ Extraction radius 

* Approximate percentage of total signal 
Total number of extracted counts. 

Estimated number of background counts. 
Minimum number of counts per spectral bin. 
Number of spectral bins. 



22 


Table 5 
Spectral Fit 


source 

number 

model" 

xl 

V 

Nh/1022 

-2 

cm 

r or kT’’ 

Radius" 

km 

fW/IO-'^ 

ergcm“^s“' 

19 (PSR) 

pi 

0.89 

31 

0 16^®'®^ 
'^■^"-0.08 

1 

^'^^-0.13 


o C4+O.O6 
'^•'^^-0.05 

26 

bb 

1.10 

44 

0.13!^^ 


1 97+0.29 
^•^^-0.23 


26 

pi 

0.86 

44 

0.68![!;“ 

c 04+0.39 


14.8!°:| 

26 

nsa 

0.96 

44 

^•^^-0.04 

'^•'^^-0.01 

14.5!|^ 

3-35!j:? 

26 

mekal 

0.88 

44 

n 072+0.02 

^••^•^-0.05 

0.59;°:“ 


3.26!}:“ 

4 

bb 

1.39 

30 

0.16![!;j« 

1.07!0;P 

0.063!n? 

1.69!°:“ 

4 

pi 

1.14 

30 

0.86!«:“ 

1 cq+0.15 
^•'^^-0.24 


9 40+0.71 
^•^^-0.26 

4 

mekal 

1.14 

30 

0.77![!;j9 

i3.3;^V 


9 94+0.14 
^'•^^-0.48 

17 

bb 

1.73 

21 

0.0-0.028 

0.88!«:«« 

0.073!ni 

1 07“'*° 

17 

pi 

1.02 

21 

0.17;0;P 

1-30!^:?^ 


1 47+0.15 

17 

mekal 

0.99 

21 

Q la+O.ll 
^-0.17 

oq+42 

•^^-26 


1 44+0.13 
^•^^-0.19 

28 

bb 

0.67 

15 

0.64![!;| 

0.67!°:P 

0.095!°:“ 

0.59!^:“ 

28 

pi 

0.74 

15 


3.08!!;“ 


2.20!°:11 

28 

mekal 

0.68 

15 

1 OO+0.63 

^■bo_Q47 



1.09!°:“ 

25 

bb 

1.43 

8 

0.0-0.13 

l-12![lf 

0-026!0:|i 

0.340!°:°' 

25 

pi 

0.77 

8 

0.0-03 

0.82!°;« 


U.^Z2)_qq2 

25 

mekal 

1.11 

8 

0.17!«i 

14.2-79.80 


u.:)/d_oo3 


Note. — 


" bb = blackbody; pi = power law; nsa = neutron star H-atmosphere; mekal = optically-thin 
thermal plasma. 

* The entry in this column depends on the spectral model — it is the power law index F or 
the temperature kT in keV {kT^ for the nsa model). 

Blackbody radius (7?^ for the nsa model). 

Unabsorbed flux in the 0.5-8.0 keV band. 



